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ABSTRACT 


This paper proposes an analysis framework for robustness analysis of a nonlinear dynamics system that can be rep- 
resented by a polynomial linear parameter varying (PLPV) system with constant bounded uncertainty. The proposed 
analysis framework contains three key tools: 1) a function substitution method which can convert a nonlinear system 
in polynomial form into a PLPV system, 2) a matrix-based linear fractional transformation (LFT) modeling approach, 
which can convert a PLPV system into an LFT system with the delta block that includes key uncertainty and schedul- 
ing parameters, 3) //-analysis, which is a well known robust analysis tool for linear systems. The proposed analysis 
framework is applied to evaluating the performance of the LPV-fault detection and isolation (FD1) filters of the closed- 
loop system of a transport aircraft in the presence of unmodeled actuator dynamics and sensor gain uncertainty. The 
robustness analysis results are compared with nonlinear time simulations. 

Nomenclature 


a 

angle of attack (AOA), (rad, deg) 

A 

elevator deflection, (rad, deg) 

Q 

pitch angle rate, (rad/s, deg/s) 

S a 

stabilizer deflection, (rad, deg) 

V 

true airspeed, (m/s) 

T 

thrust, (N) 

e 

pitch angle, (rad, deg) 

Pa 

air mass density, (kg/m 3 ) 

h 

altitude, (m) 

Psch 

scheduling parameters 

7 

flight path angle, (rad, deg) 

Psys 

system parameters 

c 

mean chord length, (to) 

q 

dynamic pressure, (A/to 2 ) 

m 

total mass, (kg) 

Si 

the i th parameter of the LPV model 

C7 

inertia coefficient, 1 / I yy , (/cg _1 ?Ti~ 2 ) 

s 

coefficient matrix 

$zeng 

summation of z positions of engines, (to) 
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Aerodynamic coefficients 


CD Mach '■ drag coefficient at fixed mach number 
Cl ■ total lift coefficient 

CL baai c : lift coefficient for the rigid airplane at zero stabilizer angle 

C m : pitch moment coefficient 

C rnbastc : pitch moment coefficient for the rigid airplane at zero stabilizer angle 

K a : effective factor of the elevator 


1 Introduction 


In general, robustness of a nonlinear system is an open problem. Due to the large diversity of nonlinear models and 
computational difficulties, there is no unified analysis method for all possible classes of nonlinear dynamics. This 
paper is concerned with robustness analysis of a class of nonlinear systems described by polynomial functions. This 
class of nonlinear systems can be represented by a polynomial linear parameter varying (PLPV) system, whose system 
matrices are polynomial functions of scheduling parameters. Conventionally, stability and performance analysis of 
an LPV system has been formulated into a linear matrix inequality (LMI) feasibility problem or a performance index 
minimization problem with LMI constraints, which represent a Lyapunov function type of analysis [1,9, 12, 13,21,28]. 
For LMI approaches, an LPV system can be described as a set of convex sets [10] or as a set of linear systems defined 
at grid points over the scheduling parameter spaces [22,28]. Robustness analysis of an LPV system, which can capture 
a nonlinear system, unmodeled dynamics, real parameter uncertainty, and time-varying uncertainty, can be included 
in an LMI optimization [1, 12, 19] with scaling matrices. In this paper, instead of using a convex set defined by vertex 
systems or a set of linear systems defined at grid points, nonlinear system variations due to scheduling parameter 
change are converted into an LFT model with a structured delta block, which contains scheduling parameters. The key 
difference between an LMI approach and the proposed analysis framework in this paper is that the nonlinear dynamics 
are converted to an LFT model with the structured delta block, which can cover all possible variations of scheduling 
parameters. The variations defined by the structured block are treated as an uncertainty set in the robustness analysis 
framework. This analysis framework is applied to the robustness analysis of the LPV-FDI filters integrated with the 
closed-loop system of a transport aircraft. 

Robustness analysis of an integrated fault tolerant control (FTC) system with fault detection and isolation (FDI) 
filters has been of interest to the NASA Aviation Safety Program (NASA AvSP) [23]. Specifically, a fault tolerant 
control law prevents fault propagation and reduces closed-loop performance degradation due to faults. It makes the 
detection of faults by FDI filters more difficult in the integrated system. In the validation and verification process of 
a fault tolerant flight system, the robust performance analysis of FDI filters is important before implementing these 
into real systems. In [1 1, 29], robust FDI filters are designed for an linear time invariant (LTI) system under feedback 
control using LMI techniques to maximize fault sensitivity performance and to simultaneously minimize the influence 
of unknown inputs, which represent dynamic effects of uncertainty. The singular value /i analysis method has been 
used to analyze robustness of FDI filters on the LTI systems in [11], In [23], robustness analysis of an integrated sys- 
tem with an FTC law and FDI filters was presented using // analysis. The nonlinear longitudinal motion of a transport 
aircraft was converted into a PLPV system using the function substitution method. A function substitution method has 
been introduced in [27] to convert a nonlinear system into a quasi-LPV model. The function substitution method has 
been modified in [8,22] to minimize variations of decision variables over the scheduling parameter space using linear 
optimization. The method has been sucessfully applied to develop an LPV model from a nonlinear system including 
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look-up table aerodynamic coefficients [22] or polynominal fitted aerodynamic coefficients [23]. For nonlinear dy- 
namics described in polynominal form, the function substitution method is automatically coded to provide a PLPV 
model [24] using the symbolic toolbox in MATLAB. In this paper, a PLPV system is constructed from the combined 
model of the LPV model of a tranpsort aircraft and its LPV-FDI filters since the scheduling parameters for the aircraft 
model and its filters are the same. 

In [23], the nonlinear dynamics of a transport aircraft were converted into a PLPV system without considering 
the LPV-FDI filters. The PLPV system was converted into an LFT model with the delta block, which contained all 
dynamic purturbations due to changes in flight condition. In [23], robustness analysis of the integrated system with 
an FTC law and LTI-FDI filters in the presence of unmodeled actuator dynamics was presented. In this paper, a ro- 
bustness analysis framework of polynomial LPV-FDI filters is introduced. Robustness analysis of the transport aircraft 
(a Boeing 747-100/200 aircraft) is performed using three key tools: 1) function substitution method, 2) LFT model- 
ing approach, and 3) p analysis tool. The previously developed robustness analysis framework in [23] is extended 
to analyze the performance of a class of nonlinear multi-variable PLPV dynamics systems. The robustness analysis 
framework is illustrated by analyzing the robust performance of the integrated LPV-FDI filters and FTC system for a 
Boeing 747-100/200 aircraft. 

The organization of the paper is as follows. The matrix-based LFT modeling approach [5,7] is discussed in Sec- 
tion 2. In Section 3, the analysis framework for a nonlinear multi-variable polynomial system is presented. In Sec- 
tion 4, the analysis method is introduced and illustrated via an example of robustness analysis of two LPV-FDI filters 
for a Boeing 747-100/200 nonlinear longitudinal dynamics model. The nonlinear longitudinal dynamics equations 
of the airplane, the actuator models, the LTI-FTC controller, and the FDI filters are from [8], [18], and [26], and are 
presented in Sections 4.1 to 4.3. The nonlinear equations of longitudinal dynamics of the airplane and of the LPV-FDI 
filters are rewritten as polynomial functions of the deviations of the states and inputs from their trim values. In Sec- 
tion 4.4, the polynomial models of the airplane and the filters are combined and converted into an LFT model with the 
parameters in the diagonal “uncertainty” block being velocity and angle of attack. The resulting LFT model is then 
connected with the controller, the linear actuator models, and the weighting functions. The closed-loop interconnected 
system and the weighting functions are presented in Section 4.5. Results of the robustness analysis and nonlinear 
simulation are presented in Section 4.6. Section 5 presents concluding remarks. 


2 Preliminary Theory Background 


There are a few methods [16, 17] of obtaining an LFT model from a multi-variable PLPV model. One of those methods 
is the matrix-based LFT modeling approach [5-7, 14]. The approach converts multi-variable PLPV models to their 
equivalent LFT models. Some of the key concepts of the matrix-based LFT modeling approach are summarized in 
this section because this approach is used in this paper, and also because these key concepts give insights into the 
equivalence between a multi-variable PLPV model and its LFT model. 

Consider 


X 



A(c>i, . . . , Si, . . . , 5„ ) B(Si, <)„) 


X 

— S(p sys ) 

X 
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C(5i, ...,Si,,..,S n ) D(Si, ...,S i: ..., S n ) 


u 


u 


where S, contains a constant nominal component <5, o and a bounded varying component (i, A 

Si = S io + Si A (2) 

Psys represents the vector of system parameters ..., Si, ..., S n , and A(p sya ), B(p sya ), C(p sya ), and D(p sya ) are 
multi-variable polynomial matrix functions of Si. It is apparent that S(p aya ) is also a multi-variable polynomial matrix 
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function of 5-i. If the coefficient matrix of the term S^ 1 . . . S'-' . . . SIR, where r/i is the exponent of Si, is denoted by 
S s vi 5 Vi sH n > S(Psys) can be expressed as 


m i rrii m n / n \ 

S(Psys) = E • • • E • • • E W 

rji—O rji—0 r] n = 0 \ i—1 / 

where m,, which is finite, is the maximum degree of <5,; . The coefficient matrix Sgo s o s o , denoted as S a , is the 
nominal component of the LPV model. 

Consider a sketch of an LFT model in Fig. 1. The matrix A(p sys ) contains the system variation, and can be 
represented as follows 

A( ps-ys ) — tliag[(liT^ 1 , . . . , Silf ti , . . . , S n Ik n ^ (4) 

The structures of P, R, L, Q matrices, suggested by [5], are 
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The LFT equation associated with Fig. 1 is given below 


R = 


Ri 


Rn 


(5) 


(6) 


S(p sys ) = S A ( Psys ) + S 0 = L(I - AP)~ 1 AR + Q = LA(I - PA^R + Q. (7) 

If (I — A P)~ 1 (or (/ — PA) -1 ) is replaced with a Taylor series expansion, the right hand side of Eq. (7) has the form 
of the right hand side of Eq. (3). With the given coefficient matrices, matrix theory is used to solve Eq. (7) for P, R, 
L, and Q such that the following conditions are satisfied: 1) P is a nilpotent matrix so that the Taylor series expansion 
of the matrix inversion is finite, and 2) p, P,j, Ri (i = 1 , ,n, j = 1 , ,ri) are found such that the expansion of 
L(I — AP) -1 AP + Q matches with the given coefficient matrices term by term. The size of the identity matrix R. 
associated with S. t does not have to be known; it is found when the equation is solved. The dimension of the system 
variation block A(p sys ) is k, . Details of how Eq. (7) is solved are beyond the scope of this paper and can be 
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found in [5-7]. Other LFT modeling methods are based on a symbolic approach [3, 17]. It is noted that the LFT model 
equivalent to the given LPV model is not unique. The system variation block of the LFT model found by any of the 
three methods does not necessarily have the smallest dimension possible. 

The parameters in A are normalized to be consistent with the convention of p analysis. The normalization of the 
parameters can be done before or after the LFT model is obtained. Before normalization, the range of parameter Si is 
[5i min , assuming Si is bounded. After normalization, the range of Si is [—1, 1]. Although it seems trivial to 

mention the fact that — 1 < normalized <5, < 1, this remark helps explain the results of /; analysis in Section 4. The 
maximum singular value of the normalized variation block (<r(A)) can be thought of as “distance” from the nominal 
of the normalized parameters; thus, the // value is the inverse of the maximum “distance” from the nominal of the 
normalized parameters, which still satisfies the performance of the system. 

In brief, this section presents some key concepts in the matrix-based LFT modeling approach to transform a given 
multi-variable PLPV model into an LFT model. A multi-variable PLPV model and its equivalent LFT model, there- 
fore, are two mathematical ways of expressing the same model. Thus, // analysis on the LFT model yields another 
approach to analyzing PLPV models. 


3 Analysis Framework 


Consider a nonlinear closed-loop system, that has the form of Eq. (1), where A , B, C, and I) are polynomial matrix 
functions of p sys . A robustness analysis problem is defined as “Can the closed-loop system achieve a pre-defined 
performance level in the presence of constant but unknown uncertainty such as unmodeled actuator dynamics and 
sensor gain uncertainty over the operating envelope?” In the analysis framework, there are three steps: 

1. The nonlinear system is converted into a PLPV system using a function substitution method [8,23,24], 

2. The PLPV system is converted into an LFT model with a structured variation block of the system scheduling 
parameters using the matrix-based LFT modeling approach described in the previous section. 

3. The augmented closed-loop system with performance weighting matrices, actuator unmodeled dynamics, and 
sensor gain uncertainty in the form of an LFT model is analyzed using the /j analysis tool. 


Robust performance analysis of LPV-FDI filters using this analysis framework are presented in the next section. 


4 Robust Performance Analysis of LPV-FDI Filters of Boeing 747-100/200 


The FDI filters of the airplane are considered robust in this paper if they can detect faults over the considered flight 
envelope of [150, 250] m/s for velocity and [—2, 8] degrees for angle of attack at the considered altitude of 7000 m. The 
variation block of the system parameters can be thought of as a structured uncertainty block, so the word “robust” here 
agrees with its conventional meaning. The robustness of the integrated LPV-FDI filters and closed-loop FTC model in 
the conventional sense (i.e., in the presence of model uncertainty) is also considered in Section 4.6. For simplicity, it is 
assumed that the airplane model and the actuator models have no uncertainty unless mentioned otherwise. Fig. 2 shows 
the block diagram of the closed-loop control of Boeing 747-100/200 and its FDI filters. In the sketch, p sc /, represents 
the scheduling parameters. The nonlinear closed-loop FTC model and the LPV-FDI filters are taken from [8, 18,26], 
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commands 



Figure 2: Block diagram of closed-loop system of a Boeing 747-100/200 aircraft and its FDI filters. 


4.1 LTI controller and actuators 


The LTI controller of the system in Fig. 2 is an controller, which is from [18]. The controller has eighteen states, 
six inputs, and three outputs. The inputs of the controller are 7 cm d — 7 , V /g, 9, q, V, V crn( i, where the subscript 
cmd stands for command and (7) denotes the deviation of a state from its trim value. The outputs of the controller are 
signals to the elevator actuators, stabilizer actuators, and the engines (thrust actuators). Following [26], the transfer 
function of the elevator is ^g 7 , and the transfer function for thrust is - + Q 5 5 , The stabilizer input to the airplane will 
not be used in the robust performance analysis of the FDI filters because the stabilizers are assumed to stay at their 
trim positions. The stabilizer transfer function is, therefore, not included. 


4.2 Nonlinear model of Boeing 747-100/200 

The nonlinear longitudinal dynamics model of Boeing 747-100/200 are from [ 8 ]. 


a = [1 


qSc 
2 mV 2 


(1.45- 


l-8a: cg ) 


dC L 


qS „ dC L 


dq ]q+[ ~mV Kl 


dS e 


]S e + [- 


4 

— — (sin a + 0.0436 cos a)]T + A 
mV 


( 8 ) 


A = 


1 

V 


(sin a sin 9 + cos a cos 9)g 



(9) 


q = CT 2V° i^tr ~ |( X - 45 ~ l-8x c .g.)^(eosax c .g. + sin a z c . g .)]q 

+c 7 qScK a [ ^ - |^f(cos ax c .g. + sin a z c . g )\5 e ( 10 ) 

+ c 7 qScK a d c™ S s + c 7 S zeng T + Q 

Q = c 7 qScC mbasio + c 7 qS[C DMach (cos a z c , g , - sina5 c . ff .) - C Lhasic { cosax c . g . + sina^ c . g .)] (11) 


• L ± 

V t = — (cos a 
m 

V = (sin a cos 9 — 


— 0.0436 sin a)T + V 
cos a sin 9)g - —Cr> Mach 


( 12 ) 

(13) 
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(14) 

(15) 


9 = 


q 


h = (cos a sin 9 


sin a cos0)V 


The altitude of the airplane is weakly coupled with other states via air mass density, and speed of sound, which affects 
the aerodynamic coefficients. It is assumed that the air mass density and speed of sound are constant in the neighbor- 
hood of the considered altitude. With these assumptions, the state h can be removed without affecting other states in 
Eqs. (8-15). It is noted that state removal is problem-dependent. 

Before the above highly nonlinear equations are approximated as polynomial functions of velocity and angle of 
attack, several sub-problems must be solved. The aerodynamic coefficients are tabulated in look-up tables [18], and 
thus must be polynomial fitted. The polynomial fitted aerodynamic coefficients, using the method in [15,20], in the 
velocity range of [150, 250] m/s, and the angle of attack range of [—2, 8] degrees at the altitude of 7000 m are 


ach 

k 20 °? + KiqCX + KqiV + ftoo 

(16) 

dC L 

d5 e 

= T02V 2 + TqiV + Too 

(17) 

C L bas . 

(c = ? 7 io a + 7701V + 7700 

(18) 

>asic 

^20 a 2 + Cioct + C01V + ^00 

(19) 

dC m 

dS e 

= C02V 2 + C01V + Coo 

(20) 


where 


«2o = 3.2653 

t 02 = -1.35 93 x 10 -7 

ryio = 5.1806 

£20 = 2.3923 

C02 = 4.3535 x 10" 7 


K10 = 3.4772 x 10 -2 
r 01 = 3.9396 x 10" 5 
7701 = 1.2560 x 10" 3 
60 = -1.4614 
Cot = -1.1644 x 10" 4 


K01 = 4.4526 x 10" 5 
too = 3.5418 x 10" 3 
Voo = -4.0826 x 10" 3 
£0-1 = -3.2036 x 10" 4 
Coo = -1-7626 x 10 -2 


koo = 9.9258 x 10" 3 


Coo = 0.1264 


( 21 ) 


Other aerodynamic coefficients are assumed to be constant [18], = 4.1, ■= —17.6, = —2.8. The 

aerodynamic coefficients are polynomial functions of V and a. The nonlinear equations Eqs. (8-15) with tabulated 
aerodynamic coefficients are called high fidelity nonlinear (HNL) equations; those equations with polynomiial fitted 
aerodynamic coefficients are called low fidelity nonlinear (LNL) equations; the LNL equations, after all the nonlinear 
terms are expressed as polynomials of the states, and the state h is removed, are called polynomial (PN) equations. 
Later, HNL, LNL, and PN will be applied to the equations of the LDI filters with the same connotation. Let x t 
be the trim value of the state x. Let x be the deviation of the state x from its trim value (x = x — x t ). The 
trim values of the states and inputs of the nonlinear model with polynomial fitted aerodynamic coefficients (LNL 
model) are [a t , qt, V, 0 t , h t \ = [1.05 deg, 0 deg/s, 227.02 m/s, 1.05 deg, 7000 m], and [ S et , S St , T t ] = 
[0.163 deg, 0.590 deg, 42291 N], Lor comparison, the trim conditions of the high fidelity nonlinear (HNL) model of 
Boeing 747-100/200 are [a t HNL, qtHNL , V t HNL, V ti hnl, 9 t HNL , h t HNi\ = [1.25 deg, 0 deg/s, 227.42 m/s, 

1.25 deg, 7000 m], and [S etHNL , S StHNL , T tHNL ] = [—3.19 deg, 1.99 deg, 42291 N], The trim conditions of the 
LNL model will be used for the analysis. The quantity 1/17 is approximated by its first order Taylor series expansion, 
1/Vt — {l/Vt)V. It is assumed that a, 9 are small angles, thus sin a ~ a, sin 9 ~ 9, cos a ~ 1, cos 9 ~ 1. After the 
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Figure 3: Time responses of the open-loop system under 
elevator input at the altitude of 7000 m. 


Figure 4: Time responses of the closed-loop system under 
7 - and velocity-commands at the altitude of 7000 m. 


polynomial approximation and function substitution, the nonlinear model, including the terms. A, Q, and V, can be 
written in the following form 

i = A(V,a)x + B(V,a)u. ( 22 ) 

The time responses of the open-loop high fidelity nonlinear (HNL), low fidelity nonlinear (LNL) and polynomial 
(PN) models under elevator input are compared in the following plots (Fig. 3). The elevator inputs are 1 degree when 
0 < t < 10 s, -1 degree when 10 < t < 20 s, 1 degree when 20 < t < 40 s and 0 when t > 40 s. The plots show 
the deviations of the states of open-loop models from the trim states of the LNL model. The time responses of LNL 
model and PN model are indistinguishable. The time responses of the closed-loop HNL, LNL and PN models under 
7 and V commands (without elevator faults and thrust faults) are compared in Fig. 4. The plots of V, 7 of the HNL, 
LNL and PN models almost overlap. Also shown in the plots are velocity command (cmd) and 7 command. Note 
that there is a small difference in the phugoid period due to the polynomial fitting of the aerodynamic coefficients 
and the approximation of the nonlinear terms. Overall, the plots show that the developed PN model is useful for later 
robustness analysis of the integrated closed-loop system. The LNL model is just a middle step to obtain the PN model. 


4.3 FDI filters 


This section briefly describes the LPV-FDI filters for the longitudinal motion of a Boeing 747-100/200 in [26], The 
elevators of the airplane are lumped together as one actuator, and all four engines are considered as another actuator, 
as assumed in [26]. Reference [26] chose to design one filter for all the elevators and another filter for all engines. The 
equations of the FDI filters in matrix form [26] are 
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( 23 ) 


tbi(f) — Ni (Psch)'Wi (f) Gi(p sc h')x(t') Fi(p sc f l ')'Uf(t') 

n{t) = MiWi{t) - Hix(t) 


Ni(Psch) ~ N.i0 + NuPsc^ + • • • + Nigp sc h a 

Gi(pach) = Gi 0 + G i\ Psch\ + ' ' ' + Gigp sc hg 
F'ii'psch) F,o 3“ F{\p sc } ll T * * * T F.igp sc h 9 


(24) 


where vj t is the state vector of filter i, x is the state vector of the airplane model, Uf consists of the desired (unfailed) 
actuator outputs (<5 e , S s , T) (see Fig. 2) and a constant input to the filters (uf = [<5 e , 5 S , T , 1]) [26]. The index * = 1,2 
refers to the FDI filters. Filter 1 is the elevator FDI filter; filter 2 is the thrust FD1 filter. Filter 1 has four states, and 
filter 2 has three states [26]. Nij, G, ? - , F , 3 , Mi, //,, i = 1,2 , j = 0, . . . , 9 are constant matrices as given in [26]. 
The notation p sc i, :i , j = 1 ..... 9 denotes scheduling parameters. The scheduling parameters are defined as p sc h 7 = q, 

q 1 n q dCj, q - dCm - 

Psch 2 — yi Psch-a — y, Pschi — u — a, p S ch 5 d ibasic V’ Pschg — gg e y > Psch 7 — '~ D rnac h Q’ Pschs — gs c 9’ 

Pschg = G mbasio q in [26], The aerodynamic coefficients are approximated as polynomial functions of V and a, so the 

scheduling parameters are functions of V, a and 6 of the airplane. 

If the nonlinear FDI filters are to be put in polynomial form, the same sub-problems as those of the airplane model 
must be solved, i.e., the polynomial fitting of the aerodynamic coefficients, the approximation of the nonlinear terms, 
and the untilization of the function substitution method. As with the airplane model, the trim values of w, must be 
obtained. At trim flight conditions without actuator faults, u = Uf (see Fig. 2), r.i(t) = 0 by the design of the FDI 
filters, the states of the FDI filters are constant, and the rates of the states of the FDI filters are zero. Let Wi t denote the 
trim values of the states of the i th filters, let Wj be the deviations from the trim values, and let f j(<) be the deviation 
of the output of filter i from its trim value ry t . The trim values of uq are found by setting x, Uf, Wi to Xt , Ut, 0, 
respectively, and solving for tv,, . Eq. (23) can be rewritten as 


m(t) = Ni{V, a)wi(t) - Gi(V,a,9)x(t ) + Fi(V,a)uf(t) 

n{t)=fi(t)+ n t = MiWi(t) - Hix(t) + MiW it - HiXt = MiWiit) - Hix(t) ^ 25 - ) 

o o 

where N-i (x) , Gi ( x ) and F i: (£■) are N t , Gi and F t after function substitution. With the assumptions of constant air mass 
density and constant speed of sound in the neighborhood of the considered altitude, it was noticed that the altitude of 
the airplane, the fourth state of the elevator FDI filter, and the third state of the thrust FDI filter were weakly coupled 
with other states of the FDI filters and had negligible effects on the output of the FDI filters. Those states are removed 
from the PN models of the FDI filters. As mentioned previously, the state removal is problem-dependent. 

The time responses of the two FDI filters under elevator fault and thrust fault are presented in Fig. 5. The legend 
“HNL” indicates that the HNL model of the airplane and HNL model of the FDI filters are in the closed-loop system. 
The legend “LNL” indicates that the LNL model of the airplane and the LNL model of the FDI filters are in the closed- 
loop system. The legend “PN” indicates that the PN models of the airplane and the FDI filters are in the closed-loop 
system. The simulation is started with no faults and zero initial states of the FDI filters. The dashed lines in Fig. 5 
represent time-varying actuator faults. Piecewise constant elevator faults are assumed to occur at 50, 80, 120, 160, 200 
and 250 seconds, and thrust faults are assumed to occur at 40, and 150 seconds. The initial conditions of the HNL FDI 
filters are not at trim values [18,26]; therefore, the time responses of the HNL filters have an initial transient. It takes 
about 20 seconds for both HNL filters to converge to the faults, and after that the time responses of all the models of 
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the FDI filters are close to each other. It is observed that for an elevator fault of 2 degrees from 120 t/l second to 1 60* / ' 
second, the output of the elevator FDI filter is about 2.6 degrees (see Fig. 5); that is, the elevator FDI filter has about 
one-degree fault estimation error per three-degree elevator fault. The observation will be used later to determine the 
weighting function for the performance of the elevator FDI filter. 
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Figure 5: Time responses of the LPV-FDI filters under elevator faults and thrust faults. 


4.4 Combining the models of Boeing 747-100/200 and the FDI filters 

Rather than considering models of the airplane and the FDI filters separately, they are augmented because they have 
the same scheduling parameters (see Eq. (22) and Eq. (25)). The combined state vector is x p = [x: w\\ W 2 ], where x, 
w\, W ‘2 are the states of the corresponding PN models. The combined equations of the plant and filters becomes 


x 


* 


A(p sys ) 

0 

0 

B(psys) 

0 


U>1 

= 

— Gl(f) S y S ) 

Nl(p s y S ) 

0 

0 

Fl(psys) 

< 

W 2 { 


. —G2{p s ys) 

0 

-^2 (/^sys) 

0 

F 2 i.fi sys) _ 



V 

s 


W\ 

W 2 

U 

Uf 




(26) 


p sys refers to V and a (9 in Eq. (25) is treated as a state variable rather than a scheduling parameter). The vectors u 
and Uf are inputs to the airplane model and inputs to the FDI filters, respectively. Inputs to the airplane are the actuator 


outputs plus faults, if any (Fig. 2). Inputs to the FDI filters are the predicted actuator outputs. S (see Appendix) is a 
polynomial matrix function of V and a. The highest degree of V in .S' is 4 and that of a is 1 . Eq. (26) is expressed in 
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LFT format using the matrix-based LFT modeling approach as in Eq. (1). The variation block is 

A(V,a) = diag[V I nv , aI Ha \ (27) 

where ny = 11 and n a = 2. The range of V(= V — V t ) is [—77.0, 23.0] m/s. The range of &(= a — at) is 
[—0.053, 0.121] rad. The LFT model is normalized based on the ranges of V and a to obtain A(V r ,ci). Like A, 
A is not the “uncertainty” block, but represents the nonlinear variations of S with respect to the normalized V and 
normalized a. 

Certainly, the LFT model of the airplane and the LFT model of the LPV-FDI filters can be obtained separately and 
connected together to form an LFT model of the airplane and filters. However, one of the advantages of LPV-FDI 
filters is that their dynamics vary with the dynamics of the airplane, in this case, via V and a. Separate LFT models 
for the airplane and the filters lead to independent variation blocks for the airplane and the filters, and this does not 
utilize the strength of the LPV-FDI filters. Comparison of the analysis results of the integrated system using separate 
LFT models and combined LFT models will be presented in Section 4.6. 



Augmented system for robustness analysis 

Figure 6: Closed-loop interconnected system and weighting functions. 


11 




4.5 Closed-loop system with weighted performance 


The closed-loop system with weighted performance is described in Fig. 6. The block “B747-100/200 Aircraft & FDI 
filters” with the feedback block “A” is the normalized LFT model in Fig. 1 . The transfer functions T d d , T ldv and 
weighting functions in Fig. 6 are 


Tid = 

diag(Tid^ , Tid v ) 

Tid P = 

0.35 2 

Tid v = 

0.15 2 

s 2 +0.7s+0.35 2 

s 2 +0.3s+0.15 2 

w p = 

diag(W Pe , W Vv ) 

w Pe = 

^0(^ + 1) 

w pv = 

lOCKifo+l) 

0.005 + 1 

0.005 “l" 1 

w pr = 


W pr 1 = 

0.33(0. ls+1) 2 
(2s+l) 2 

II 

<N 

S- 

0. 001(0. 001s+l) : 
(ios+ 1) 2 


W f = diag { 3 deg, 10000 N) 


T t d 7 and T idv , taken from [26], are the desired tracking models from 7 cm d to 7 and from V cm d to V, respectively. W Pe 
is the weighting function to penalize the flight path angle tracking error. W Pv is the weighting function to penalize the 
velocity tracking error. W pr 1 and W pr 2 are the weighting functions on the difference between the elevator fault and 
the output of filter 1 and on the difference between thrust fault and the output of filter 2, respectively. W pr 1 allows 
about one -degree error in the detection of three-degree elevator faults at low frequency. Recall that the output of the 
elevator FDI filter is about 2.6 degrees for a two-degree elevator fault (Fig. 5). Therefore, for a three -degree elevator 
fault, the error in the fault estimate is about 1 degree. W pr 2 allows 1000 N error in the detection of 10000 N thrust fault 
at low frequency. At high frequency, W pr 1 and W pr 2 relax the weights on the differences. W c is given in Section 4.6. 


4.6 Results of robust performance analysis 


This section presents the robust performance analysis results of the integrated LPV-FDI filters and the closed-loop 
FTC system of a Boeing 747-100/200 at 7000 m over the considered flight envelope, V = [150, 250] m/s and 
a = [—2, 8] deg. This section also compares the performance analysis by using the combined LFT model and the 
separate LFT models of the airplane and the LPV-FDI filters. In addition, the performance of the linearized LTI-FDI 
filters at the LNL trim conditions is considered. Nonlinear time simulations are performed on the HNL models to 
validate the /./, analysis results. In this analysis, the weighting matrix block W c (see Fig. 6) describes the maximum 
size of the allowable commanding signals. Three levels of commanding signals are considered; they are (case 1) 
W c = diag( 1 deg, 1 m/s), (case 2) W c = diag(3 deg, 5 m/s), and (case 3) W c = diag ( 6 deg, 10 m/s). 

The weighted nominal performance from the four inputs, 7 command, velocity command, elevator fault and thrust 
fault, to the four outputs, 7, velocity, elevator residual and throttle residual, is presented in Fig. 7. The word “nominal 
performance” is borrowed from the robust control jargon; it means the variation block A in Fig. 6 is zeroed out (V = 0 
and a = 0), and A„ and A s are zeroed out. The H a c norm must be less than or equal to one at any frequency. If 
Hoq > 1 at some frequency, the nominal system does not satisfy the defined performance criteria. In that case, further 
robust performance analysis is not necessary. The plot in Fig. 7 shows the nominal system performs well at the LNL 
trim conditions for (case 1) commanded V signal of 1 m/s, commanded 7 signal of 1 deg, (case 2) commanded V 
signal of 5 m/s, commanded 7 signal of 3 deg, and (case 3) commanded V signal of 10 m/s, commanded 7 signal 
of 6 deg. It is observed from Fig. 7 that the nominal performance gets worse as the command signals increase. The 
performance of the integrated system under varying command signals was studied. Other performance analyses in this 
section focus on the performance of the integrated system under the same bounded command signals, but with varying 
uncertainty models. 

Because the H ^ norm is less than one in the considered frequency range, the performance of the integrated LPV- 
FDI filters and the closed-loop FTC system at any combinations of V and a is considered. The normalized variation 
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Figure 8 : Performance of the system over the flight enve- 
Figure 7: Nominal performance: case 1 (dash dot); case 2 lope without uncertainty: case 1 (dash dot), case 2 (dash), 
(dash); case 3 (solid). case 3 (solid). 


block A is non-zero, but A u and A s are still zeroed out. The inputs of the system, 7 command, V command, elevator 
fault and thrust fault, and the outputs of the system, W Pe , W Pv , W pr 1 , and W pr 2 , are connected by fictitious full 
uncertainty blocks, and the plots of the upper bound for /r are presented in Fig. 8 for three cases of the command 
gains. The word “performance” in the titles of the sub-plots in Fig. 8 refers to the outputs of W p in Fig. 6 . The top left 
subplot shows the performance of the airplane without faults, the top right subplot shows the coupling effect from the 
V and 7 commands to the outputs of the FDI filters, the bottom left subplot shows the coupling effect from the faults 
to the airplane performance, and the bottom right subplot shows the performance of the FDI filters with zero command 
signals. From the top left subplot of Fig. 8 , /i > 1 for case 3 reveals that the airplane cannot satisfy the performance 
requirements over the considered flight envelope when the commanded V signal is 10 m/s, and commanded 7 signal 
is 6 deg. Even though the airplane performance degrades, the performance of the LPV-FDI filters is satisfactory over 
the considered flight envelope under the same command signals as shown in the right subplots of Fig. 8 . 

This method does not only reveal how the system performs at trim conditions or how the system performs if the 
parameters are within their ranges, but it also helps answer many questions regarding the performance of the LPV- 
FDI filters. For example, the same method can be used to study the performance of the LPV-FDI filters over new 
ranges of the considered parameters ( V and a for this particular problem), or how the filters perform with a slightly 
different model of the airplane, or how the filters perform when the system model has some uncertainty. Several of 
the above questions are considered here. It is assumed that the uncertainty of the actuator models are represented as 
post multiplicative uncertainty. Let W urlf and W unt denote uncertainty models of the elevator actuator and throttle 
actuator, respectively. 

The bounds for /j for case 1, shown in Fig. 9, are calculated in the presence of one percent V and a sensor gain 
uncertainty (dash-dotted), 1 . 2 % (dashed) and ten percent (solid) unmodeled actuator dynamics over the considered 
flight envelope. Fig. 9 shows the performance of the airplane and the LPV-FDI filters under the influence of the 
command signals and faults. The top right subplot and bottom left subplot of Fig. 9 reveal that there are coupling 
effects from the command signals to the filter residuals and from the faults to the airplane performance. Fig. 10 shows 
the performance of the LPV-FDI filters under the combined effect of the command signals and faults. Note that the 
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1 . 2 % uncertainty actuator models are selected as 


^un e — ^urt f 


0.012(s/0.5 + l) 
s /100 + 1 


(29) 


The gain 0.012 is selected to make the upper bound for /i for the performance of the FD1 filters close to one but not 
more than one. Eq. 29 represents the maximum unmodeled actutator dynamics that still guarantees satisfactory filter 
performance for case 1 over the flight envelope (/x < 1, see the dashed lines in Fig. 9). The dash-dotted lines (/< < 1) 
shows that the integrated system performs well for case 1 in the presence of one percent sensor gain uncertainty. The 
solid lines in Fig. 9 show that the LPV-FDI filters cannot achieve the pre-defined performance requirements even for 
case 1. The lower bound for /j in Fig. 10 confirms that // > 1 for the integrated system with ten percent unmodeled 
actuator dynamics. This means that there exists a flight condition within the considered flight envelope, where the 
performance of the LPV-FDI filters degrades. In Fig. 10, the values of // are less than one for one percent sensor 
gain uncertainty, close to one for 1 . 2 % unmodeled actuator dynamics, and greater than one for ten percent unmodeled 
actuator dynamics. 

The HNL dynamics of the integrated system are simulated with the LPV-FDI filters and the FTC law for selected 
command and fault signals (shown in Fig. 1 1 and 12) to validate the robustness analysis results. For case 3 without 
uncertainty (A u = 0, A s = 0), velocity is commanded as a 10 m/s increase at 50 s and a 10 m/s decrease at 180 
s with a 10 m/s velocity command offset at 0 s. The velocity command offset leads the aircraft to flight condition 
of V = 238 m/s, which is close to the boundary of the considered flight envelope. It is noted that the output of the 
elevator LPV-FDI filter is bounded within 1 deg, and the output of the thrust LPV-FDI filter is bounded within 1000 
N as allowed by the weighting functions in Eq. (28). The time responses of the EDI filters (marked with the symbol 
* in Fig. 1 1 and Fig. 12) show that the LPV-FDI filters achieve the pre-defined performance requirements at the flight 
condition close to the boundary under commanded V signals of 10 m/s and commanded 7 signals of 6 deg. The 
satisfactory performance of the LPV-FDI filters even for case 3 confirms the robustness analysis results in Fig. 8 . The 
lines with the symbol y show the time responses of the integrated system with 1 . 2 % unmodeled actuator dynamics. 
The commanded 7 signal is 1 deg up at 20 s and 1 deg down at 100 s, and the commanded V signal is 1 m/s increase 
at 50 s and 1 m/s decrease at 180 s with velocity command offset of 10 m/s. The lines with the symbol A show the 
time response of the system with the same pattern of the command signals in the presence of one percent sensor gain 
uncertainty with the same velocity command offset. The outputs of the LPV-FDI filters, corresponding to /j, < 1 cases, 
stay within bounds. The time responses of the system agrees with the robustness analysis results in Fig. 9, which 
predicts a satisfactory performance of the FDI filters in the presence of one percent sensor gain uncertainty or 1 .2% 
unmodeled actuator dynamics. The lines with the symbol • in Fig. 11 and Fig. 12 show the time responses of the 
integrated system with 7 command 1 deg up at 20 s and 1 deg down at 100 s and velocity command 1 m/s up at 50 s 
and 1 m/s down at 180 s with zero offset. In this case, the LPV-FDI filters do not satisfy the performance requirements 
since the time responses of the residual signals are out of bounds. Again, the nonlinear time simulation correlates to 
the robustness analysis results in Fig. 9. Note that all vertex points of the uncertainty blocks A,, and A s are considered 
for the nonlinear simulation for /r < 1 cases; for // > 1 cases, only one time history of the filter outputs at one vertex 
point is shown (see Fig. 1 1). 

The performance of the FDI filters was considered at different flight conditions, under different command signals, 
for different values of /i, and for different types and different levels of uncertainty. However, it is of interest to 
investigate how the FDI filters perform under the combined effect of sensor gain uncertainty and unmodeled actuator 
dynamics. The performance of the system with one percent sensor gain and one percent unmodeled actuator dynamics 
is shown in Fig. 13 and Fig. 14. Fig. 14 indicates that the FDI filter performs well because fi < 1. The outputs of 
the FDI filters, shown in Fig. 15, at the flight conditions, presented in Fig. 16, stay within bounds. Again, the time 
simulation correlates to the // analysis results. 

For comparison of the separate LFT models, the ft bounds of the system using the separate LFT models are 
calculated and shown in Fig. 17. Note that the LPV-FDI filters allow the dynamics of the filters to vary with the 
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Figure 9: The upper bound for // of the system with one 
percent uncertainty on the gains of the velocity sensor and 
the angle of attack sensor (dash dot), with 1 .2% uncertainty 
on the actuators (dash) and with ten percent uncertainty on 
the actuators (solid) for case 1. 


Figure 10: The upper bound for /i of the system with one 
percent uncertainty on the gains of the velocity sensor and 
the angle of attack sensor (dash dot), with 1.2% uncertainty 
on the actuators (dash), the lower and upper (solid) bounds 
for [i for the system with ten percent uncertainty on the 
actuators for case 1 . 



Time (s) 






Figure 1 1 : Time history of the LPV-FD1 filters for the sys- 
tem with no uncertainty for case 3 (*), with one percent 
uncertainty on V and a sensor gains for case 1 (A), with 
1.2% uncertainty on the actuators for case 1 (y), and with 
ten percent uncertainty on the actuators for case 1 (•). 


Figure 12: Time history of the closed-loop HNL system 
for the system with no uncertainty for case 3 (*), with one 
percent uncertainty on V and a sensor gains for case 1 (A), 
with 1 .2% uncertainty on the actuators for case 1 (y ), and 
with ten percent uncertainty on the actuators for case 1 (•). 
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Faults to performance 



Frequency (rad/s) 


Commands to residuals 
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Frequency (rad/s) 
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Figure 13: System with one percent uncertainty on the 
gains of the velocity sensor and the angle of attack sen- 
sor and one percent uncertainty on the actuators dynamics 
for case 1 . 


Figure 14: Performance of the FDI filters with one percent 
uncertainty on the gains of the velocity sensor and the angle 
of attack sensor and one percent uncertainty on the actua- 
tors dynamics for case 1 . 




Time (s) 




Time (s) 
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Figure 15: Time history of the LPV-FD1 filters for the sys- 
tem with one percent uncertainty on V and a sensor gains 
and one percent uncertainty on the actuators for case 1. 


Figure 16: Time history of the closed-loop HNL system 
with one percent uncertainty on V and a sensor gains and 
one percent uncertainty on the actuators for case 1. 


16 


dynamics of the model of the airplane. That is the reason why the model of the airplane and the LPV-FDI filters are 
combined into a plant before the LFT model is computed for the combined plant. The parameters in the variation 
block of the LFT model insure that the LPV-FDI filters still capture the dynamics of the airplane after the combined 
plant is put in LFT form. If the LFT models of the airplane and the LPV-FDI filters are found separately and then 
connected without combining the parameters in the variation blocks, the parameters in the variation block of one LFT 
model are treated by the // analysis as independent from the parameters in the variation block of the other LFT models. 
That does not utilize the strength of the LPV-FDI filters because the dynamics of the LPV-FDI filters do not depend 
on the dynamics of the airplane. Also, the // analysis result of the system for case 1 is more conservative, as shown 
in the right subplots in Fig. 17 because the parameters in each variation block are independent from parameters in 
the other variation blocks. The variation blocks of the components such as the aircraft dynamics, the LPV-FDI filter 
(filter 1) for elevator residual, and the LPV-FDI filter (filter 2) for thrust residual are A a / C = diag( V ' / 7 , ci/ 2 ), 

A f uteri = diag(Vl 5 , a), and A futer 2 = diag(Vl 2 , a), respectively. For /./, analysis, eachblockis combined into 

A = diag( A a / C , A f ut. eri , A futer 2 )- The combined block contains 14 copies of V and 4 copies of a. Recall that 

the A block in Eq. (27) contains 1 1 copies of V and 2 copies of a for the LFT model of the combined system with 
the LPV-FDI filters and the aircraft model. It is observed from the left subplots in Fig. 17 that the upper bounds for // 
shown in the left subplots are equal to the // values shown in the left subplots in Fig. 8 since the // values in the left 
subplots are not related with the FDI filter dynamics. 

For comparison of the linearized FDI filters, the lower and upper bounds for \i of the system with the linearized 
FDI filters are calculated for case 1 and are shown in Fig. 18. It is observed from the right subplots in Fig. 18 that the 
linearized FDI filters do not satisfy the pre-defined performance criteria over the flight envelope. Fig. 8 and Fig. 18 
show that the LPV-FDI filters perform well in a bigger flight envelope than that covered by the linearized LTI-FDI 
filters. 

Evaluating the performance of FDI filters with nonideal sensors (sensors with noise) and subjected to disturbance 
is a topic for future study. Another topic for future investigation is reconfiguration of the control laws using FDI filter 
outputs to make them robust in the presence of system failures. Extension of this analysis to other control methods is 
also of high interest. 


5 Conclusions 


The robustness analysis framework is presented for polynomial nonlinear systems using three key tools: 1) function 
substitution method, 2) LFT modeling approach and 3) /t analysis tool. The key of the analysis method is the LFT 
representation of a nonlinear multi-variable polynomial system, which can capture the nonlinear dynamic variation 
due to flight condition change. In the framework, it is easy to integrate with unmodeled actuator dynamics and sensor 
gain uncertainty. The analysis framework is applied to the robustness analysis of the integrated LPV-FDI filters and the 
closed-loop FTC system of a Boeing 747-100/200. Robustness analysis results are validated via nonlinear simulation 
in the presence of unmodeled dynamics and sensor gain uncertainty. 
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Appendix: The coefficient matrix S 


The polynomial coefficient matrix S with all the states of the airplane model and LPV-FDI filter models and all the 
control inputs are included below. Deleting the rows and columns corresponding with the altitude, the fourth state of 
the elevator FDI filter, the third state of the thrust filter, and the stabilizer, shifting the last three columns from the sixth 
row to the last row to the right, and adding zeros where appropriate yields the same S in Eq. (26). 
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